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Abstract 



Qh! We investigate the universal low-energy dynamics of the simplest Higgs 

I phase for gravity, 'ghost condensation.' We show that the nonlinear dynam- 

ics of the 'ghostone' field dominate for all interesting gravitational sources. 
I Away from caustic singularities, the dynamics is equivalent to the irrota- 

^ I tional flow of a perfect fluid with equation of state p oc p^, where the fluid 

particles can have negative mass. We argue that this theory is free from 
catastrophic instabilities due to growing modes, even though the null energy 
condition is violated. Numerical simulations show that solutions generally 
have singularities in which negative energy regions shrink to zero size. We 
exhibit partial UV completions of the theory in which these singularities are 
smoothly resolved, so this does not signal any inconsistency in the effective 
theory. We also consider the bounds on the symmetry breaking scale M 
in this theory. We argue that the nonlinear dynamics cuts off the Jeans 
instability of the linear theory, and allows M < 100 GeV. 
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1 Introduction 



Is general relativity the correct description of gravity at long distances and times? 
Certainly, there are good reasons for thinking that this is the case. Experimentally, 
gravity has been probed at distance scales ranging from 10~^ mm (in short-range 
force experiments) to at least 10^^ cm (the size of the solar system). Theoretically, 
general relativity is the unique Lorentz-invariant theory of massless spin 2, and its 
conceptual elegance is beyond question. However, the situation is far less clear on 
cosmological distance and time scales. Structure formation, galaxy rotation curves 
and gravitational lensing, and the accelerating expansion of the universe cannot be 
explained by general relativity coupled to known matter. These anomalous effects 
are conventionally attributed to 'dark matter' and 'dark energy.' However, given the 
fact that the observed effects are purely gravitational, it makes sense to ask whether 
they may have a common origin in a modification of gravity in the infrared. These 
considerations have led to a revival of interest in consistent infrared modifications of 



In the present paper, we further investigate the model of Ref. [5], 'ghost conden- 
sation.' This can be viewed as the universal low-energy dynamics associated with the 
simplest Higgs phase for gravity, arising when Lorentz symmetry is broken sponta- 
neously. Breaking of Lorentz symmetry is of course ubiquitous. For example, time- 
dependent fields in cosmology define a preferred frame. However, these solutions are 
not the ground state of the theory: they carry energy density and dilute away as the 
universe expands. Any 'modification of gravity' induced by such solutions becomes 
relevant only at scales of order the horizon. Wc are instead interested in the situation 
where Lorentz symmetry is broken in fiat spacetime, allowing nontrivial modification 
of gravity inside the horizon. This means that the symmetry breaking sector has 
peculiar properties; in particular, the stress-energy tensor must vanish in the ground 
state: 



Spontaneous breaking of Lorentz symmetry gives rise to a gaplcss scalar excitation 
analogous to the Nambu-Goldstonc bosons that arise from spontaneous breaking of 
internal symmetries. 'Ghost condensation' gives rise to a single such mode, and is in 
this sense the minimal model of spontaneous breaking of Lorentz symmetry. We refer 
to the scalar mode as a 'ghostone boson.' In analogy with the Higgs phase for gauge 
theory, the ghostone mode mixes with the graviton, modifying gravity in a nontrivial 
manner. 



gravity [1, 2, 3, 4, 



5, 6, 7]. 



(1.1) 
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Ref. [5] studied this theory and analyzed the modification of gravity in the weak- 
field fimit. The dynamics is governed by a consistent effective theory defined by the 
scale M where the symmetry is broken. It was shown that the ghostone mode gives 
a possible new origin for dark energy and dark matter. Ref. [8] showed that the 
ghostone mode may also be the infiaton, leading to interesting testable consequences. 

In the present paper, we further investigate the dynamics of this theory. We show 
that nonlinearities dominate the dynamics of the ghostone sector for all gravitational 
sources of interest. In particular, the time scale for the onset of nonlinear dynamics 
for a fixed gravitational source is precisely the infall (or orbit) time associated with 
the source. Away from singularities, the nonlinear solutions are equivalent to the 
gradient flow of a fluid with equation of state 



The gradient fluid flow picture breaks down at caustic singularities where geodesies 
of fluid particles cross. At these points higher-derivative terms become important. 
Numerical simulations show that these remove the caustic singularity in most cases, 
but in some cases singularities still form. These correspond to regions of negative 
energy that shrink in size until their dynamics becomes sensitive to the underlying 
microscopic theory. We argue that these can be stabilized in a partial UV completion 
of the model into 'gauged ghost condensation.' 

We also address stability issues in the nonlinear theory. Although the energy is not 
positive-definite in the nonlinear theory, there are no growing modes corresponding 
to infinite positive energy radiated to infinity. We argue that the worst instability is 
that of negative energy regions shrinking in size, as found in the simulations. 

Based on these results, we address the question of the bounds on the ghostone 
sector. We find that the most sensitive bounds come from the fact that the Jeans 
instability found in Ref. [5] produces regions of positive and negative energy and 
that those regions induce scintillation of light rays coming from far distances by 
gravitational lensing. Demanding that the total angular deviation by the scintillation, 
or the random walk, not exceed the observed angular resolution gives the bound 
M < 100 GeV. Other possible observable effects such as modification of gravity near 
astrophysical sources and energy loss of moving sources are argued to be negligibly 
small. Black holes in the ghost condensate are studied in Ref. [9], where it is found 
that they also do not give strong constraints on the model (see also Ref. [10]). 

This paper is organized as follows. In section 2, we review the Higgs phase for 
gravity, emphasizing the universal nature of the low-energy dynamics. We also review 
the origin of the relevant length and time scales in the linearized theory. In section 
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3, we discuss the nonlinear dynamics of the ghost condensate. We discuss the time 
scales and present the fluid picture. In section 4, we show the numerical simulations of 
the nonlinear evolutions and discuss the resolution of caustic singularities. In section 
5, we discuss the bounds of this theory. This includes mass and energy accretion 
in slow-moving objects, gravitational lensing and energy loss from moving objects. 
Wc do not claim to have a complete understanding of the dynamics, so this section 
is intended to be preliminary and provocative. In section 6, we briefly discuss the 
possibility that ghost condensate may constitute the dark matter. We show that the 
initial growth of the density perturbations in the linear regime is identical to that of 
the standard cold dark matter. Whether it can form the correct structure depends 
on the details of the nonlinear evolution which is left for future investigations. Our 
conclusions are presented in section 7. 

2 Review of the Linecir Theory 

2.1 Effective Theory 

What is a Higgs phase for gravity? It is easiest to answer this question in linearized 
general relativity, where we expand the metric about flat space 

Qu-u = Vix,^ + V (2-1) 

and keep only terms quadratic in /i^^. In this theory, the fields hf^u are closely analo- 
gous to gauge fields with gauge transformation law 

Sh^.^-{d^C. + duQ, (2.2) 

where are the generators of infinitesmal diffeomorphisms 

x''^x^ + C". (2.3) 

We want to consider the case where the time diffeomorphisms generated by are 
spontaneously broken. This means that time diffeomorphisms are realized nonlinearly 
in the effective theory containing only the ghostone field. The minimal model contains 
a single real ghostone field tt that shifts under time diffeomorphisms: 

Stt = -Co. (2.4) 

Note that tt naturally has units of time. We now write the most general effective 
Lagrangian invariant under these symmetries. This contains the Einstein Lagrangian, 
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and additional terms constructed from the invariants 



E = TT — ^hoQ, 

Kij = 7:ihij - 9ihoj - djhoi + 2c?j5,-7r). 



(2.5) 



(2.6) 



The leading terms are 




(Note that we have built in the fact that fiat space is a solution by not writing any 
linear terms in the Lagrangian.) 

In the limit where we turn off gravity, we see that the ghostone mode has dispersion 
relation 



where a = cti + 0:2. If the dimensionless couplings in the effective Lagrangian are 
order 1, the mass scale M is the scale of new physics in this theory. 

Eq. (2.7) is the starting point for analyzing the linear dynamics of the theory, 
including the coupling to gravity. However, there are several points that may not be 
completely clear in this formulation. First, it may not be completely clear in what 
sense the breaking of Lorentz symmetry is spontaneous, since the time direction 
appears to be explicitly singled out. Second, it is not clear whether the absence of a 
linear term in S (allowed by all symmetries) requires an additional fine-tuning, like 
the cosmological constant which is a linear term proportional to /i^. Third, we would 
like to have an efficient way to extend this to the full nonlinear theory. All of these 
issues can be elegantly understood in the covariant formulation, which we turn to 
next. 



Consider an effective theory with a real scalar (j) that is invariant under a global shift 
symmetry 



where A is a spacetime independent constant. For example, (f) could be a Nambu- 
Goldstone boson from the breaking of a global symmetry, or a 0-form gauge field in 
string theory. 



2.2 Covariant Formulation 



54) ^X, 



(2.9) 
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A conventional effective Lagrangian for tliis tlieory is 



C = (2.10) 

witli equation of motion 00 = 0. Tliis lias solutions with = n^jX^ for any constant 
4- vector n^. If is timelike, we can choose the time direction so that the solution is 

= (2.11) 

This is a solution for any constant c. At first sight it may appear that these are 
obviously not candidate ground states of the theory, but the situation is actually 
more subtle. Suppose we expand in fluctuations about this solution 

= ct + 7r. (2.12) 

Note that under time diffeomorphisms, tt transforms as 

^TT = -Co/c, (2.13) 

just like the ghostone mode considered above. Expanding the Lagrangian to quadratic 
order in tt, one finds that the fluctuations for tt have good time and space kinetic 
terms. This means that the solution is stable under local fluctuations for any value 
of c ! The reason for this is that the theory has a conserved current associated with 
the shift symmetry 

= 9^0. (2.14) 

Solutions with c ^ have a constant nonzero charge density. Local excitations 
cannot change the total charge, so configurations with lower energy cannot be reached. 
However, when we turn on gravity, solutions with c 7^ will cause the universe to 
expand, and the charge will dilute away. Lorentz invariance is therefore not broken 
spontaneously in this theory. 

Consider instead an effective Lagrangian of the form 

£efr = M4p(X), X^d^<\>d^<\>. (2.15) 

Note that has dimensions of length (or time), so that X is dimensionless. This 
omits only terms with more than one derivative acting on 0, such as (D^)^. We 
will include their effects below. This also has solutions of the form = for any c. 
Expanding in fluctuations Eq. (2.12), we obtain 

Aff = { \lc^P'\c^) + P\c^)\ e - P'{c'){V7iy + C(7r3)} (2.16) 
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We see that small perturbations are stable provided 

2(PP"{(P) + P'{^) > 0, P'(c2) > 0. (2.17) 

(Note that a conventional kinetic term P{X) = +|-^ satisfies these conditions.) 

To see whether any of these solutions may be regarded as candidate ground states, 
we compute the stress-energy tensor: 

T;, = [-P{X)g^, + 2P'(X)9>9,0] . (2.18) 

In the solution = ct, the first term is proportional to g^i, and can be cancelled by 
tuning the cosmological constant. The second term gives rise to 

roo = 2c2p'(c2), Tij^O. (2.19) 

This gives rise to an expanding universe unless 

P'{(^) = 0. (2.20) 
Note that the conserved shift current is given by 

= 2P'{X)d^<P, (2.21) 

so the charge density vanishes in this configuration. Therefore, if P{X) has the form 
shown in Fig. 1, then the configuration with P'{X) = is a candidate ground state. 
This is the 'ghost condensate.' From now on, we rescale the field 4> so that the ground 
state is at X = 1, and write 

E = i(X - 1). (2.22) 

Stability in the linearized theory requires higher-derivative terms to give a nonzero 
spatial kinetic term, for example 

A£efr = —{n<py = ^(W)' + ---. (2.23) 

This leads to the o;^ ~ k'^ dispersion relation for the ghostone mode already discussed 
above. 



6 



P(X) 




X 



Fig. 1. The kinetic fuction P{X) for ghost condensation. 



2.3 Higgs Phase for Gravity 

The ghostone mode mixes with gravity, giving rise to a nontrivial modification of 
gravity in the infrared. A detailed analysis is given in Ref. [5], but much of the 
physics can be understood from very general considerations, which we review here. 

Recall first the ordinary Higgs mechanism for gauge theory. In this case, a nonzero 
charge gives rise to a tadpole for the Goldstone mode. Excitations of the Gold- 
stone mode themselves carry charge, and because opposite charges attract, the Gold- 
stone charge tends to screen the original charge. If the charge configuration is time- 
dependent, the Goldstone field responds at the speed of fight. This is illustrated 
schematically in Fig. 2. 

Similarly, in the Higgs mechanism for gravity, a nonzero mass (the charge that 
sources gravity) gives rise to a tadpole for the ghostone mode. Excitations of the 
ghostone mode themselves carry energy, but because positive masses attract in grav- 
ity, the ghostone field an^z-screens the source gravitational field. This gives rise to 
an instability of the vacuum analogous to the Jeans instability of pressureless mat- 
ter coupled to gravity. The ghostone Jeans instability can be seen in the dispersion 
relation for the scalar mode that follows from the mixing with gravity: 

This is illustrated schematically in Fig. 3. The length and time scale associated with 
the Jeans instability are 



Mpi Ml 
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Gauge Higgs Phase 




Screening 

Exponentially falling 
Potential 



Fig. 2. A schematic illustration of the screening effect for gauge theories in the Higgs 
phase. 

For M < 10 MeV, Tj is longer than the lifetime of the universe and there is clearly 
no constraint from the Jeans instability. 



2.4 Negative Energy 

Even for time shorter than the Jeans time scale T7, stability is an issue because the 
ghostone energy can be negative. Recall that the ground state X = 1 is the boundary 
of the stability region. Expanding to higher orders in tt, we find interaction terms 
such as 



-'eflf 



(2.26) 



In regions where tt < the cubic term gives rise to negative gradient energy, and 
therefore an instability. Classically, the (V^tt)^ term restores stabihty at short wave- 
lengths, and there is a growing mode only for long wavelengths: 

A > (2.27) 

Quantum-mcchanically, there are fluctuations at all length scales. Nonetheless, Ref. [5] 
showed that there is no quantum instability in the effective theory using a scaling ar- 
gument. If we scale energy hj E ^ sE, the quadratic kinetic terms are left invariant 
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Gravity Higgs Phase 




Fig. 3. A schematic illustration of the anti-screening effect for gravity in the Higgs 
phase 



if we scale 



t s'H, 

-1/2. 



(2.28) 



With this scaling the cubic operator 7r(V7r)^ in Eq. (2.26) scales as s^^'^ and is therefore 
(barely!) irrelevant. All other operators arc even more irrelevant, showing that there 
is a regime of low energies where the expansion is under control. 



2. 5 Preview of Nonlinear Effects 

The arguments above show that the effects of the nonlinear terms are under control 
at low energies and small field amplitudes. However, in the presence of large classical 
gravitational sources the nonlinear terms can become important. In fact, the time 
scale for the ghostone field near a classical source is just the gravitational infall time 
of the source. This can be understood from the form of the stress-energy tensor. In 
the approximation where the Lagrangian is £ = Pi^), the stress-energy tensor has 
the form 

oc -PiX)g^, + 2P'{X)u^u,. (2.29) 
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where 



= 9^0. (2.30) 

This has the form of a stress-energy tensor for a perfect fluid with 4- velocity m^.^ 
The fact that is a gradient means that the flow is irrotational: di^Ui,] = implies 

— * _ 

Vxu — 0. The equations of motion for the ghostone field follow from the conservation 
of the stress-energy tensor V^T^^, — 0, which are intepreted as the conservation of 
energy and momentum in the fluid. In the presence of a classical gravitational source, 
a fluid clearly responds on a time scale given by the infall time, and therefore so does 
the ghostone mode. As it will be shown in the next section, this is exactly the time 
scale where the nonlinear term becomes important. 

For small fluctuations about the minimum X — 1, P{^) can be approximated by 

p{x)^\{x-if = h:\ (2.31) 

where an overall contribution to the cosmological constant has been removed. We 
can read off the equation of state of the fluid from the energy-momentum tensor. It 
is 

P=4rr (2.32) 

where 

p = M^E. (2.33) 

For p <^ M^, small ghostone fluctuations therefore behave like an almost pressureless 
fluid, which makes the ghostone excitations a candidate for the dark matter in our 
universe. We will have more to say on this later. 

As we noted above, Hubble expansion drives E — > 0, and therefore p — > 0. Even 
for p = the fluid has nontrivial dynamics in the presence of gravitational sources: 
the fluid particles follow geodesies. As a result, the local preferred frame (where 

= d^(f) is 'at rest') is a freely-falling frame. This frame is also the preferred frame 
from the spontaneous breaking of Lorentz invariance, so we conclude that the 'aether' 
is dragged by the local gravitational field. 

Since geodesies generically cross (especially in the presence of gravitational sources), 
the field configuration develops caustic singularities where the field gradients diverge. 



^Thc defining property of a perfect fluid is that at each point there is a frame in which the 
stress-energy tensor has the form T^i, = diag(— 
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Such phenomena were also observed in some other scalar field theories [11]. Near 
these singularities, the higher-derivative Q;(V^7r)^ term can no longer be neglected. 
This term contributes positive gradient energy, and we find in numerical simulations 
that it generically smoothly resolves the caustic singularity, giving rise to a 'bounce' 
with outgoing tt waves and positive and negative regions of E near the would-be 
caustic. 

The fact that S (and hence p) can be negative brings up again the question of 
the stability of the theory. As discussed above, regions with E < modes with 
sufficiently long wavelengths are unstable. In numerical simulations, we find that 
negative energy regions tend to shrink while the amplitude of E grows inside the 
region. This can be understood from the fiuid picture, since this is valid in the limit 
where we neglect the (V^tt)^ term, which is a good approximation away from caustic 
singularities. In this picture the E < region consists of fluid particles with negative 
mass. It is therefore clear that energy (mass) cannot flow across the boundary of 
the E < region, since the boundary consists of particles with vanishing mass. The 
boundary can move however, and in the E < region the positive pressure favors 
large gradients and causes the E < region to shrink. Numerical simulations show 
that some E < regions continue to shrink until they exit the regime of validity of 
the effective theory. These singulartics need to be resolved in a more fundamental 
theory. Similar conclusion was also obtained in Ref. [12] which studied the two- 
dimensional case. However, we show that the total energy inside the E < regions 
formed in astrophysical situations is very small, and does not lead to any observable 
consequences provided the singularities are regulated in a smooth way. We will discuss 
a partial UV completion to do this, and present numerical evidence that it works. 

The nonlinear dynamics affects the bounds on M derived in the linear theory. 
The Jeans instability in the linear theory gives a bound M < 10 MeV if we require 
that there is no exponential growth of the oscillatory potential within the age of 
the universe. However, the nonlinear dynamics is expected to cut off the instability, 
and may weaken this bound. We argue below that the strongest bound comes from 
gravitational lensing due to regions of positive and negative energy produced by the 
Jeans instability. Demanding that the random walk of light rays due to the lensing 
does not smear out the observed CMB anisotropies gives the bound M < 100 GeV. 

We also consider other possible bounds on the ghost condensate from the grav- 
itational sector. The modiflcations of the gravitational potential are small due to 
velocity effects [13, 14]. We also consider energy loss in the nonlinear theory, as 
well as energy stored in would-be caustic singularities. Wc find that these effects are 
negligible, and we believe that the theory is safe for M < 100 GeV. 
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3 Nonlinear Dynamics 



We now turn to the nonlinear dynamics of the theory. The nonlinear dynamics is very 
rich and complex, and we emphasize that we do not claim a complete understanding 
in this work. It is therefore important to keep in mind that there is a simple limit of 
this theory, independent of the details of the nonlinear dynamics [5]. The Ghostone 
sector naturally couples to matter only through gravity. (Gravitationally induced 
direct couplings to standard model fields are easily seen to be negligibly small.) The 
maximum value of the gravitational energy in the Ghostone energy is of order M*^, 
which does not affect even cosmology if M < (Mpii^o)^''^ ~ 10~^ eV. Such low values 
of M are still very interesting for cosmology, since the ghost may be a source of dark 
energy and dark matter [5] and may drive inflation [8]. 

Another general point to keep in mind in the following is that the modifications 
of gravity vanish in regions where T, — {X — 1) and we neglect the (V^tt)^ term. 
This is because in this limit, the Lagrangian in 'unitary gauge' (f) — t {(f) — 0) is 

C^^gP{g'y (3.1) 

This has the form of a gauge- fixing term for the gauge X = g^^ = 1, so if X = 1 
initially, the gravitational dynamics is unchanged. This means that any modification 
of gravity will be suppressed by the small amplitude E or by the Q;(V^7r)^ term which 
is small at long wavelengths. We will see that this suppresses many possible effects 
compared to naive expectations. 

In the remainder of the paper, we consider M S> 10~^ eV to see what modifications 
of gravity may be observable today, and determine the experimental limit on the scale 
M. 

3.1 Effective Lagrangian 

We are interested in the gravitational dynamics of large, slow-moving sources, where 
relativistic and retardation effects are neghgible. We further restrict ourselves to 
the domain of weak gravity, so wc arc not considering black holes. (Black holes in 
the presence of the ghost condensate are discussed in Ref. [9].) It is convenient to 
incorporate these approximations from the beginning in an effective Lagrangian. The 
fact that gravity is weak means that we use the linearized approximation for gravity. 
The fact that the system is nonrelativistic means that dt -C V. For the motions of 
objects, we can neglect gravitational radiation, so we can parameterize the metric by 
the scalar potential $ as 

ds^ = (1 + 2^)df^ - (1 - 2$)df ^ (3.2) 
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where 

'PI 



^"^-7^- (3-3) 



We also assume that ghostone amphtudes are smaU in units of M. Since only deriva- 
tives of TT are meaningful, this means that 

IVttI < 1. (3.4) 

With these approximations, the deviation of X from its minimum is given by 

E= 1) =7r-|(V7r)2-$. (3.5) 

Note that there is nothing in the assumptions above that determines the relative size 
of the terms on the right-hand side. The effective Lagrangian is then 

/:eff = M^{lE^-^(VV}- (3.6) 

All the terms omitted are smaller than the ones we have kept by the assumptions 
above. Note that all the nonlinear effects in this approximation are due to the (Vtt)^ 
term in E. 



3.2 Nonlinear Time Scale 

We now estimate the time scale for the nonlinear gradient terms in Eq. (3.5) to 
dominate the dynamics. We assume that the gravitational potential $ is determined 
by external sources, and varies on a length scale L in the region of interest. (For 
example, at a distance r outside a gravitational source, L ~ r.) If initially tt = 0, the 
gravitational potential gives a tadpole that forces vr to be nonzero. We want to know 
the time scale Tnl when the nonlinear term in E becomes as important as the linear 
TT term. This is determined by 



71 TT^ 



T^-^, (3.7) 



where tt is a typical ghostone amplitude. This gives 

Tk. ~ A. (3.8) 

We recognize this as the gravitational response time associated with the potential $. 
For example, for a point source $ ~ i?s/?" where Rs is the Schwarzschild radius, and 
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we have T^-^ ~ r^^'^/Rs, which is the Kepler relation. This is a very direct way of 
seeing that the nonlinear effects become important on the gravitational time scale. 

Solving Eq. (3.7) for tt, we obtain 

IVttI ~ y ~ < 1. (3.9) 

That is, ghostone amplitudes are small for weak gravity, in agreement with our as- 
sumptions. 

On the other hand, in the linear approximation the Q;(V^7r)^ becomes comparable 
to TT^ at TLin ~ ML^. As a result, the nonlinear evolutions completely dominates for 



Tnl ^ Trir,. ~ ML\ or 



Z (3.10) 

For example, outside a gravitational source of mass Mgrc, this condition is 

1 /Mpi^ ^ 



and the earth's surface gravity is in the nonhnear regime for M ^ 10~^ eV. 

As we have emphasized, the modification of gravity is suppressed by E or by the 
Q;(V^7r)^ term. However, we will argue below that with the nonlinear evolution, these 
also become important at the time scale Tnl found above because this is the timescale 
for the formation of caustic singularities. This is to be compared with the time scale 
for the modification of gravity in the linear regime, which is the time scale for the 
Jeans instability. If wc take M < 10 MeV (as in Ref. [5]), the Jeans time scale is 
longer than the age of the universe, and the gravitational time scale of any object is 
a much shorter time scale. We see that the nonlinear effects completely dominate the 
dynamics in all regimes of interest. 



3.3 Fluid Picture 

We now show that the nonlinear equations of motion have a simple interpretation 
in terms of a perfect fluid if we neglect the Q;(V^7r)^ term in the action. As already 
mentioned in subsection 2.5, we can see this connection directly in terms of the stress- 
energy tensor for the ghostone mode. Near X — 1 we can approximate P{X) by 

P(X) = |(X - 1)2 = iE^. (3.12) 

Neglecting the ^(V^Tr)^ term, the stress-energy tensor has the form 

T^, = {-iS^^^, + Eu^u,} , (3.13) 



14 



where 



= 9^0. (3.14) 

Note that is nonzero and timchke everywhere. This has the form of the stress- 
energy tensor for a perfect fluid with 4- velocity u^. Because is the gradient of 
a scalar, the flow of the fluid is irrotational. The conservation of the stess-energy 
tensor V^T^i, — gives the equation of motion for the ghostone field, and also gives 
the Euler equation for the fiuid. For E <^ 1 we can read off the density and pressure 

p^M% P-'2M'^"-^4- (3-15) 

This establishes the equivalence of the ghostone theory without the Q;(V^7r)^ term to 
the irrotational flow of a perfect fluid ^. 

It is also insightful to understand the equivalence directly in terms of the equa- 
tions of motion. For the ghostone held, the equations of motion have the form of a 
conservation law 

E = V-[EV7r], (3.16) 

where S is the charge density. For a fluid made of particles carrying the conserved 
charge, the current is J = Hv, so we identify 

v^-Vn. (3.17) 

Again we see that the fluid flow is irrotational. In this fluid picture, the equations of 
motion are satisfied simply due to the fact that the fiuid particles carry their charges 
with them. It remains only to satisfy the relation between E and tt: 

E = 7r-i(V7r)2-$. (3.18) 

Taking the gradient of both sides and using Eq. (3.18) gives 

^ = _V($ + S), (3.19) 
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The Jeans length Lj = 2-iTCs/^'iTTGNP corresponding to this equation of state is Lj ~ Mpi/M^, 
where Cg = ^/dp/dp = pjM^ is the sound velocity. Intriguingly, this agrees with the Jeans length 
(2.25) in the linear theory up to a constant of order unity. This equation of state ignores the A;^ 
term, while the linear analysis does not take into account the nonlinear term Y? in p. Hence, it 
is not a priori clear whether these two Jeans length should be the same or not. Nonetheless, they 
agree. 
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where 



(3.20) 



is the time derivative along the particle worldline (also called the convective or Lagran- 
gian derivative). Eq. (3.19) is just Newton's law for a particle moving in a potential 
$ + S. Using the identifications Eq. (3.15) we can write 

-VE = --Vp, (3.21) 
P 

which shows that Eq. (3.19) is precisely Euler's equation for the fluid. ^ 

The equivalence between the ghostone dynamics and the dynamics of a perfect 
fluid was derived here in the nonrelativistic approximation and for linearized gravity. 
In fact, as shown in the appendix, it holds in full nonlinear general relativity. 

We can now understand better the justification for neglecting the a(V^7r)^ term, 
which was used to derive the fiuid picture. As discussed in the introduction, outside 
a gravitational source, the fiuid (and therefore the ghostone field) responds on a time 
scale of order the gravitational infall time of the source. For reasonable values of 
M, this is much shorter than the time for the Q;(V^7r)^ term to become important. 
However, the particle trajectories in a perfect fluid tend to cross, leading to caustic 
singularities. (For E 7^ 0, the pressure may stop these singularities, but they certainly 
occur for E = 0.) At the caustic singularity, derivatives of the velocity fleld blow up, 
and higher- derivative terms such as the o;(V^7r)^ term cannot be neglected. Therefore, 
the fluid picture is valid away from small regions where caustics form. 

Wc can also understand better the extent to which the locally preferred frame 
(the 'aether') is dragged by the gravitational fleld. We expect the universe to consist 
of domains where S is negligbly small and varies smoothy, separated by caustic 
regions. Away from the caustics, the locally preferred frame is a freely-falling frame, 
and in this sense the aether is dragged by the local gravitational fleld. 

3.4 Negative Energy 

We now re-examine the question of negative energy in the theory in the presence of 
the nonlinearities discussed above. Because the time translation invariance and 
shift invariance are spontaneously broken down to a diagonal time shift symmetry, 
there are different possible deflnitions of energy in this theory. Gravity couples to the 

^Note that the equivalence between the ghostone and fluid pictures is a kind of duality, since it 
exchanges a constraint equation with an equation of motion. 
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Noether charge of the original time translation symmetry, which we refer to as the 
'gravitational energy.' This is not the same as the Noether charge associated with 
the time translation symmetry that is unbroken in the vacuum, which we call the 
'inertial energy.' 

The inertial energy e is the energy associated with the unbroken time translation 
symmetry of the Lagrangian Eq. (3.6). (It is also the Hamiltonian density of the 
system.) Conservation of inertial energy states that 



£ = -V-p, (3.22) 



where 



e = + iE(V7r)^ + ^(V^^} (3-23) 

is the inertial energy density and 

= j-ETrVTT + [(V'7r)V7r - t:V{V^ti)\ | (3.24) 



P 

is the momentum density. In the linearized approximation, 

e = M- 



However, in the nonlinear theory the inertial energy is not positive definite due to 
the second term in Eq. (3.23). The energy density can be negative only in regions 
where E < 0. In fact it is easy to see that the energy is unbounded from below. For 
example, for tt = c|x| we have e = — |M'^r/. 

We might expect this theory to be unstable because any region can lower its 
energy by emitting positive energy radiation to infinity. (Because inertial energy is 
positive in the linear theory, radiation must have positive energy.) The reason this 
does not happen is that there is another conserved quantity in the theory, namely the 
Noether charge associated with the shift symmetry.^ Because tt shifts under the 
shift symmetry, the conservation of shift charge is identical to the equation of motion 
for tt: 

S = -V • J, (3.26) 
where E is the shift charge density, and the shift current is 

j^_EV7r + ^V(VV). (3.27) 



^The gravitational energy is a linear combination of the inertial energy and the shift charge, so 
the conservation of the stress-energy tensor does not give an independent conservation law. 
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If we neglect the a(V^7r)^ term, this conservation law is taken into account very 
directly in the fluid picture, to which we now turn. In this approximation, negative 
energy regions correspond precisely to regions where E < 0, so we consider such a 
region surrounded by E > 0.^ In the fluid picture, the conserved charge is carried 
by the individual particles, so the E < region consists of particles with negative 
charge, while the boundary of the E < region consists of particles of vanishing 
charge. Therefore, there can be no flux of charge across the E = boundary, and the 
total charge inside the region does not change. 

On the other hand, the S = boundary can move. Since VS points outward at 
the boundary, the equation of motion for fluid particles Eq. (3.19) implies that the 
particles on the boundary experience an inward force due to the pressure. Therefore, 
the E < region tends to shrink. 

These arguments show that the total shift charge integrated over a E < region 

Q = / rf^xE (3.28) 

does not change with time. This can also be seen from the fact that the shift current 
J vanishes on the E = boundary. The shift charge is not the same as the total 
inertial energy 

I Sxe. (3.29) 
7s<o 

However, if we neglect the a(V^7r)^ term, the flux of inertial energy across the E = 
boundary also vanishes, and therefore E also does not change with time. 

As with the fluid picture, these results hold beyond the approximations made 
here. This is discussed in the appendix. 

3.5 Caustic Solutions 

The fluid picture can be used to understand the structure of the caustic singularties 
that occur when we neglect the Q;(V^7r)^ term. We restrict attention to the case 
E = 0, where there is no pressure and the fluid particles follow geodesies. In this 
case, it is clear that there are caustics without the tt(V^7r)^ term. It is possible that 
the pressure resolves the caustics in important situations such as inside galaxy halos 
made of ghostone dark matter, but we will not consider that here. 

^Wc do not consider configurations that have E < at infinity. These can be eliminated by 
physical initial conditions. For example, in the expanding universe Hubble friction drives can drive 
S — > from above. 
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To introduce the subject of caustics, we consider the gravitational potential due 
to a uniform sphere of matter with density pQ. Inside the sphere, the gravitational 
potential is 



i.e. a harmonic potential. If we consider an initial fluid configuration where all par- 
ticles are at rest, then all fluid particles reach the center of the sphere at the same 
time! This kind of 'perfect caustic' requires very special initial conditions and sym- 
metry. Nonetheless, we will see that an analog of the perfect caustic appears to be 
responsible for the singularities that remain in the full theory even after the Q;(V^7r)^ 
term is included. 

A more realistic situation will have less symmetry in the gravitational potential 
and the initial conditions. Nonetheless, in the presence of gravitational sources the 
fluid particle trajectories tend to cross, resulting in caustic singularities. The generic 
situation is that families of trajectories cross on a dimension-2 surface. Wc want to 
zoom in on the behavior near the caustic surface, so wc can treat the problem as 
1-dimensional, and approximate the caustic surface by a plane. We also expect the 
potential to be irrelevant on the small scales where the caustic structure is important. 
This reduces the problem to a very simple 1-dimensional one. 

In the absence of a potential, the particle trajectories are just straight lines: 

X ^ Xq + Vo{xo)t. (3.31) 

Here we are labeling each trajectory by its initial position xq ai, t — 0. The initial 
velocity function vo{xo) determines the geodesies. The problem now is to flnd the 
velocity of the particle at a given point {x, t), which is a simple kinematic problem. 

To get a caustic near a; = 0, we are interested in the case where the particle 
trajectories that start near xq — are focussed toward each other: 

^;o(0) = 0, v'oiO) < 0. (3.32) 

For small Xq we can therefore approximate Vo{xo) by a linear function of Xq: 

voixo) = (3.33) 

where T is a constant that tells us how the initial velocity varies away from xq — 0. 
Solving Eq. (3.31) for xq we obtain 
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Fig. 4. The 1-dimensional 'perfect caustic' The shaded regions shows the region 
where we can expand the solution perturbatively about the perfect caustic. 



Since the velocity is constant along any trajectory, we have 

X 



v{x) = Vq{xq) 



T-t 



We can relate this to the field tt using dxi^ = —v, to obtain 



7t{x, t) 



X 



2{t-T)' 



(3.35) 



(3.36) 



Note that this solution is scale-free, since T just gives the time to the caustic. Shifting 
t ^ t — T, we obtain 



X 



(3.37) 



To understand this solution, note that the general condition for a caustic is 

dx 



dxo 



= 



for fixed t. This gives the time to the caustic for a given value of Xq as 

1 



tc(xo) = 7 



(3.38) 



(3.39) 



For vq given by Eq. (3.33) above, we see that the caustic forms at t = T for all initial 
points Xq. We call this the 'perfect caustic,' and it is illustrated in Fig. 4. 
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In a realistic case, the function vo{xo) will be more complicated. We want to 
expand about the point where the caustic first forms, i.e. the minimum of tc{xo), 
which we take to be Xq — 0. Expanded about this point, Vq takes the form 

voixo) = + + 0{L-% (3.40) 

where we have performed a boost so that fo(0) = 0. Because of the focussing of the 
geodesies, this expansion is valid only for 

t-T 



\x\ < L 



T 



(3.41) 



This region of validity is shown in Fig. 4. 

We can now follow the same steps as for the perfect caustic, keeping higher order 
terms in Eq. (3.40). We find 

xo^-T— ' , +C>(L-^). (3.42) 

and therefore 

v{x) ^ + ^ , ^ + 0{L-^). (3.43) 

As above, this is related to the ghostone field solution by d^i^ — —v. 

This result can be used to evaluate the effect of the Q;(V^7r)^ term in the action 
near the caustic singularity. In the equation of motion for tt, it contributes a term 
proportional to ad^n. This can be thought of as a restoring force, since it tends to 
smooth out the large gradients near the caustic region. Evaluating this in the caustic 
solution Eq. (3.43) gives 

= -dlv = -^j^, + 0{L-% (3.44) 

Note that a nonzero contribution to 9^7r arises only at C(L^^). In particular, there is 
no restoring force at all in the perfect caustic solution. We beheve that this explains 
the singular behavior we find in numerical simulations below. Since the d^i: vanishes 
in the perfect caustic, it is an exact solution to the full nonlinear equations, including 
the a(V^7r)^ term. This solution therefore becomes singular and exits the regime 
of validity of the effective theory. Even though the perfect caustic is a very special 
solution, it does describe the generic caustic close to the singularity. Also, the fact 
that the perfect caustic is a scale-free solution may mean that it is an attractor in 
the full nonlinear dynamics. 
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We can use Eq. (3.44) to estimate the time and distance scale where the Q;(V^7r)^ 
term becomes important. This will happen when 

«^ ^ ^- (3-45) 

For a solution that varies on a time scale T and distance scale L, the scaling in 
subsection 3.2 shows that tt ~ L'^/T^. We therefore find that the Q;(V^7r)^ term 
becomes important for 

A«S(^] , (3.46) 

where At = T — t is the time to the caustic. The distance scale where the Q;(V^7r)^ 
term becomes important is therefore 

Note that Ax, At > M'^ as long as L, T > M'^ and LjT < 1 {i.e. the system is 
nonrelativistic) , so this is within the regime of validity of the effective theory. 



4 Numerical Simulations 



In this section we describe various numerical simulations of the ghost field which allow 
us to understand the rather exotic features of its non-linear evolution. 

The test cases we first present assume symmetry to reduce the dynamics to a 
one dimensional problem, and will see that the various symmetries have different be- 
haviours. We will discuss how caustics do indeed form in the theory. As mentioned 
earlier, this is clear for a = 0, but one would naively expect non-zero a to amelio- 
rate this problem. However, as we will show, this is not the case, and for certain 
symmetries the 'perfect' caustic remains an attractor. 

As one might expect, the singular behaviour becomes less strong as one moves 
from planar through axisymmetry to spherical symmetry. Indeed without any grav- 
itational potential we will see the planar reduction exhibits singularities, while the 
spherical theory without potential does not. However, once a gravitationally at- 
tractive potential is added, all three symmetries become singular under evolution of 
regular initial data. 

Clearly assuming symmetry can lead to unphysical behaviours, whilst the phy- 
sical situation we are ultimately interested, namely structure formation, is expected 
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to have very little symmetry. Hence this section will conclude with a study of a 3-d 
numerical evolution where no continuous symmetry is present and a moving gravi- 
tational potential seeds the ghost field growth. As expected from the 1-d examples, 
we again find singular caustics do develop, and interestingly appear to take a planar 
form. 

4-1 One dimensional evolutions: Planar, axial and spherical symmetry 

We will initially consider the case of evolution of the ghost field in the absence of 
gravitational sources, seeded instead from a local perturbation in the ghost field 
itself. Then reducing to planar symmetry tt = 7r{t, r) and we may write the equation 
for the ghost decoupled from gravity in a manifestly fiux conservative form. 



where H = drir is the gradient of the ghost field, L = ^JajM ^ the second line is simply 
the definition of S and x' are the time and space derivatives of x. We then use a 
Crank-Nicholson method to evolve our initial data, which we take to be a Gaussian 
profile in tt 



with S = initially. As discussed previously, with L = E would remain zero, but 
the higher derivative term sources E. 

In order to make contact with the epoch of structure formation we wish to have 
moderate initial amplitudes so | ttq |< 1 but is still of order unity. Prom (4.2) we 
have chosen units to have initial data with unit spatial variation. We then wish to 
have L, the ghost length scale, to be L ^ 1. Naturally numerical methods limit 
our ability to separate the initial data length scale from the ghost length whilst 
maintaining accuracy. However we may separate the scale sufficiently to see the 
asymptotic properties of taking L very small. 

The first figures we show, 5 and 6, illustrate the generic behaviour for the planar 
system for small L, showing both H and S for the evolution. We show both the 
evolution of data for L = 0.005 and also for L = 0. for comparison, and the initial 
amplitude was taken to be ttq = 0.1. 

As expected L = 0. evolution leads to a caustic, whose time for formation scales as 
^caustic — 1/^0- Adding the higher derivative term changes the evolution dramatically 
around this time scale, completely smoothing out the ghost field. However, we see 



(4.1) 



7r(t = 0) = TTge 
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from the figure that the radiated waves scattered by the action of this higher derivative 
term in fact are attracted back to the origin where they grow and become singular. 
Obviously since the symmetry is planar, this growth is not due to a measure factor, but 
rather is due a 'perfect caustic' forming, which as mentioned before cannot be rescued 
by the higher derivative term we use here. Indeed this perfect caustic formation can 
be seen locally in detail from E near the singularity. 

We note that taking sufficiently large L ~ sJWq the field is completely scattered and 
the behaviour is essentially that of the linear theory and the non-linear terms never 
contribute to the evolution. However this large ghost length is clearly not compatible 
with our expectations for structure formation, where we expect L <^ \'Kq\ < 1. For 
such parameter values the behaviour shown in these figures 5 and 6 is generic, and in 
figure 7 we plot the (inverse of) the time to the perfect caustic formation compared 
to TTo for various L. As expected, this time increases with decreasing amplitude 
(eventually crossing over to the smooth linear evolution regime at L ~ when 
T — > oo), or with decreasing L (so that the radiation from the increasingly perfect 
central caustic region is slower). 

Thus for small L a cartoon of the evolution is that the initially collapsing ghost 
field causes the higher derivative term to radiate strongly, but the process of radiation 
simply refines the collapsing region into the perfect caustic form where it eventually 
collapses. 

Suitable modification of the conservative equations of motion (4.1) introduce the 
geometric measure factor associated with axisymmetry. In this case we find the 
behaviour essentially analogous to the planar case. Whilst the singular behaviour 
appears 'weaker', taking longer to reach the singularity for the same parameters, 
the singularity does indeed form. However moving to spherical symmetry we find a 
change in behaviour. 

For spherical symmetry we again take the initial data (4.2) and evolve this in the 
absence of any gravitational potential. In figure 8 we plot the evolution of H for 
TTo = 0.1 and L = 0.005 (as shown earlier for the planar case). The behaviour is 
clearly different with a totally non-singular evolution for all times, a portion of which 
is shown in the figure. Whilst for L = 0. obviously the spherically symmetric ghost 
field exhibits the usual caustic singularity, we see the higher derivative term, aided 
by the geometric measure factor, can now radiate sufficient energy to avoid any later 
energy build up. 

So far we have discussed the case of evolution of a local perturbation in the ghost 
field. We have seen that the symmetry of the initial data has a strong effect on 
whether the evolution is singular or not. Physically, however, we are interested more 
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Fig. 5. Figure on tlie left sliows H = dr-Tr for planar evolution of initial Gaussian in 
TT with ttq = 0.1 and L = 0.. Note a caustic forms as expected. The figure on the 
right shows the same evolution but with non-zero L = 0.005. Both evolutions are 
shown up to the time where the field becomes singular. We see the higher derivative 
term smooths the initial singularity present in the L = 0. case but again results in a 
later singularity. 



in the growth of the ghost field in response to a gravitational seed. Thus we must 
understand whether singularities in the ghost field form under these circumstances. 

For planar symmetry we may add a gravitational potential to the system 4.1 by 
simply introducing a potential $ to modify S as, 

H = dr S + + $ (4.3) 

We then use initial data for the evolution where both H and S are zero, and allow 
the gravitational potential to seed the growth in the ghost field. We find that for 
the weak potentials expected in our physical context, in the planar case this does 
not influence the behaviour appreciably. The gravitational source simply serves to 
produce a local non-zero configuration of the ghost field, which, once formed, then 
undergoes the singular evolutions described earlier, and is then largely independent 
of the presence of the potential. The same is true for axisymmetry. 

The interesting case is clearly spherical symmetry, as without a potential we find 
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Fig. 6. Figure sliowing S for tlie same evolution witli L = 0.005 as in figure 5. 



a non-singular evolution for smooth localized initial data. We take a potential such 
that, 

= $06^" (4.4) 

In analogy with our earlier simulations $o now acts like ttq and we have taken units 
so the physical size of the source is of order unity. Then in figure 9 we show a typical 
evolution of H and E for parameters $0 = 0.2 with L = 0.02. 

We now see a complicated evolution where waves are radiated early on as the 
higher derivative term becomes important near the symmetry axis, but later these 
waves, consisting of positive and negative S regions, are drawn grudgingly into the 
centre, eventually forming a singularity. This contrasts with the case without poten- 
tial where such waves form, but are quickly radiated to infinity. 

Thus the less singular behaviour present in spherical symmetry appears to be 
overcome once a gravitational potential is introduced and we conclude that, with the 
choice of higher derivative term in (4.1), the evolution of the ghost field in response 
to gravitational sources for all 3 reduced symmetries is singular. 



26 



0. 12 
0.1 
0. 08 
0. OG 
. 04 
0. 02 

0.02 0.04 0.06 0.08 0.1 

Fig. 7. Figure showing the inverse of the time to the 'perfect caustic' singularity, 
1/T, against initial amplitude for various L in planar symmetry. 

4-2 Three dimensional evolution 

We now test our claim that evolution is singular for gravitational seeds. In particular 
it is not obvious which symmetry will dominate the singularity evolution. For a point 
seed one might assume spherical symmetry to be the most relevant. However we have 
seen it is in the planar case that the singularity formation is quickest. 

We evolve the vr field itself in three dimensional Cartesian coordinates using a 
simple explicit second order method. Clearly here resolution becomes an important 
issue and we are less able to separate the 2 length scales in the problem. However even 
with very modest resources and resolutions (~ 150^) we can achieve useful results. 

Again we use initial conditions where vr and S are zero. We seed the dynamics of 
the ghost field by taking a Gaussian potential, 

$ = -$06""' (4.5) 

where now = (x — Xo(t)Y + {y — yo(t)Y + and, 

Xo(t) = 0.7 cost yoit) = 0.7 sin t (4.6) 

so our gravitational seed executes a circular motion in the z = plane. Hence our 
system now has no continuous symmetries although there is a preferred plane, and also 
the geometry of the source is locally spherical. Therefore we may gain information 
about the local geometry of any singularities that form. 

Figures 10, 11 and 12, show vr in the a; = plane at various times for a typical 
evolution with parameters $o = 2. and L = 0.2 (note the physical scales are not 
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Fig. 8. Plot oi H = dfTi for spherical symmetry with no gravitational potential. 
Initial data is as in the previous planar case, and again the simulation shown has 
TTo = 0.1 and L = 0.005. We see no sign of singularity formation here, and indeed 
continued time evolution confirms this. 

so well separated as in the one dimensional examples). The y = slices look very 
similar. We see that, as indicated from our study of the reduced dimension examples 
above, a singularity does indeed form, the last time frame being immediately before 
the configuration becomes singular. For the same static source, so xo(t) = yo(t) = 0, 
the singularity occurs at t ~ 3.1, but for the moving seed described here, it occurs at 
t ~ 6.2, and the source executes approximately one revolution. 

We see from the figures, and also the complete 3-d data, that the singularity 
develops in the z = plane in which the source moves, and appears to have a planar 
form, filling the disc bounding the motion of the source. This is despite the fact that 
the source only completes one 'orbit' and thus is sufficiently slow moving that this 
disc symmetry is not obviously imposed. Thus we take this as evidence confirming 
our naive expectation that the local geometry of the singularity will be the symmetry 
where singularity formation is strongest, namely planar. 



28 



Fig. 9. Plot of H = dr-Ti and S for spherical symmetry with a gravitational potential. 
We see that a singularity does eventually form, whereas earlier we saw that with no 
potential, smooth localized initial data evolves in a non-singular manner for spherical 
symmetry. 

4-3 Interpretation of numerics 

In subsection 3.4 we have shown that a negative S region does not grow but shrinks 
asymptotically. To be precise, the acceleration of the boundary of a negative energy 
region is towards the negative region itself. It has also been argued in subsection 3.5 
that the perfect caustic is an attractor. With numerical results at hand, it is easy 
to confirm them. In this subsection, for simplicity we illustrate the confirmation in 
the planar case only, while the same qualitative features are seen also in the spherical 
case with the potential. 

It is certainly worth while asking how caustics will form during structure forma- 
tion. Here, we define a caustic in the broad sense as a region where S becomes 0(1). 
Our understanding of caustics is still primitive and this question is obviously beyond 
the scope of this paper, but we can at least say that a caustic region does not swallow 
the whole space. The essential reason for this is the existence of a conserved charge 
J Tjdv resulting from the shift symmetry: if |E| grows then the volume shrinks. In 
particular, we have shown and numerically confirmed that a negative S region shrinks. 
In this sense caustics in the ghost condensate are somewhat similar to Shockwaves 
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Fig. 10. Figure sliowing vr in the x = plane at time t = 2.0 for a full 3-d evolution. 
The ghost field is sourced by a slowly moving gravitational potential, orbiting in the 
z = plane. The intersection points of the orbit with the slice shown are indicated 
by red dots in the contour plot. The peak in the field reflects the recent passage of 
the source through the plane. 

in hydrodynamics. For a shockwave in hydrodynamics, we do not need to specify 
microscopic physics, such as atomic physics, to describe its dynamics and influence 
to macroscopic physics outside. Therefore, it is expected that we should be able to 
predict what happens outside the caustic region without too many assumptions about 
a UV completion in the caustics. In the next subsection we shall discuss possible UV 
completions. 

Based on our numerical result, we expect that a planar caustic, being pretty much 
like a pancake or a disk, forms during structure formation. Indeed, the 3-d simulation 
showed that the initial stage looks spherical but deviations from spherical symmetry 
grow and the system evolves towards a planar shape. 

As future work, we should ask at what scale caustics forms first in realistic sit- 
uations. In particular, it is interesting to ask whether (i) a large caustic forms first 
and cascades to smaller ones via fragmentation or (ii) many small caustics form and 
buildup to large caustics. In either case, we might be tempted to speculate that the 
distribution of caustics should somehow trace the matter distribution since the for- 
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Fig. 11. Figure showing an intermediate time slice of tt at t = 4.0 for the same 3-d 
evolution as in the previous figure. 

mation of caustics is enhanced by the gravitational potential due to matter. Clearly, 
we need more detailed analysis to see if this is the case or not. For this purpose we 
of course need to specify the initial condition of vr and the matter distribution, which 
determines the external gravitational potential. To be specific, let us suppose ghost 
inflation [8] happened at if ^ 10~^M. This is sufficient to explain the amplitude of 
primordial density perturbations 5p/ p ~ 10~^, which determines the initial conditions 
for matter, or equivalently the gravitational potential. We also have primordial per- 
turbations of 7r. Thus, the problem to be solved is well-defined at least in principle. 
The nonlinear dynamics of the matter-vr coupled system with this initial condition 
certainly deserves further investigation. 

4-4 Possible UV completion/ extension: gauged ghost and/or new terms 

We have seen that the nonlinear dynamics leads to formation of caustics. For the 
perfect caustics, the V^vr term in the equation of motion vanishes and, thus, does not 
act as a leading spatial derivative term. This is the very reason why the perfect caustic 
does not bounce. At the same time, this means that the spatial derivative expansion 
in this background starts differently from that in the trivial vr = background and. 
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Fig. 12. Continuing from the previous figures this plot shows the time slice t = 6.1, 
just before the field becomes singular. Note the geometry of the almost singular field 
appears planar, and is roughly confined in the disc bounded by the source motion. 

thus, new terms should be included in the leading Lagrangian for the perfect caustic 
background. It is also possible to consider new degrees of freedom which become 
relevant as the system approaches the caustic singularity. 

With this in mind, for practical purposes (such as numerical studies of the ghost 
condensate as a dark matter candidate) it is probably useful to have simple models 
of new terms and/or new degrees of freedom which manifestly cutoff the singular 
behavior. Implementations adopted in this subsection are not particularly based on a 
consistent derivative expansion nor a more fundamental theory. Nonetheless, models 
here may be useful for practical purposes and, at least, show the existence of regular 
systems without caustics. 

First, we consider the introduction of new degrees of freedom by turning on the 
gauge coupling for the gauged condensation proposed in ref. [15]. (Inclusion of new 
terms will be considered later in this section.) The gauged ghost condensation has 
new degrees of freedom, but they decouple from vr in the limit of zero gauge coupling. 
Hence, by promoting the gauge coupling constant to a function of vr's derivatives it 
is possible to make the new degree relevant only near the singularity. 
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The leading Lagrangian for the gauged ghost condensation is obtained by gauging 
the shift symmetry and is of the foUowing form. 



C=-{a,-M^-—^) --F^F'^.-iv.Af. (4.7) 

In the — > hmit, F^j, goes to zero and, thus, is written as a derivative of a 
scalar: = d^n. Therefore, in this limit the original (ungauged) ghost condensation 
is recovered. In rcf. [15] the nonlinear dynamics of the gauged ghost condensation 
was analyzed. It was shown that gauging the ghost condensation resolves the caus- 
tics for the spherically symmetric case with an external gravitational force. On the 
other hand, the caustics in the perfectly planar symmetric case remains unresolved. 
It was argued that small fluctuations on top of the perfectly planar symmetric layer 
should grow and lead to fragmentation of the layer into small pieces and caustics with 
codimension more than one should not form after that. Thus, it was argued that the 
gauged ghost condensation does not lead to caustics except for extremely fine-tuned 
situations. In the following we shall implement the gauged ghost condensation as 
a possible UV completion. As a very simple implementation, the gauge coupling is 
turned on when and only when A? is large. We shall see that this indeed cuts off the 
caustics in the spherically symmetric case with external gravitational force. On the 
other hand, we do not expect to cutoff the caustics in the perfectly planar symmetric 
case since the gauged ghost condensation itself does not bounce the perfectly planar 
symmetric caustics. However, as argued in ref. [15] and stated above, a planar sym- 
metric caustics with perturbation is expected to fragment into small pieces and to 
bounce eventually. 

With the spherical symmetry, the equation of motion is reduced to 




(4.8) 



where 



S = -(X-1), e = Drd), f = ^ (4.9) 

a. 

, G = gMLphys. (4.10) 



We can eliminate E from the set of equations essentially because the electric field / 
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is determined solely by the charge. Actually, since 



= 0, 



(4.11) 



the combination in the squared bracket is independent of t. Since the expansion of 
the universe in the early universe dilutes E and / to zero with an extremely high 
precision, we set 



Hence, we obtain 



E|t=o — f\t=o — 0. 



(4.12) 
(4.13) 



and E can be eliminated from the set of equations. The result is 

1 



e = G'bf - b 



(4.14) 



We have described the gauged ghost condensation proposed in ref. [15]. Now we 
promote from a constant to a function of so that the gauge coupling is turned 
on when and only when the gradient becomes sufficiently large. For simplicity we 
choose 



Gl 



1 + tanh 
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(4.15) 



where Gq, ^0 ^-nd 5 are constants. Numerical simulation with a 



2.0 X 10-^ Gl 



1.0 X 10^, .^0 = 1-0 X 10 ^ and 5 = 5.0 x 10 is shown in Fig. 13. Fig. 14 compares 
this case with the singular result for the ungauged ghost condensation, i.e. Go = 0. 

Now, let us turn our attention to possible new terms which manifestly cutoff the 
singular behavior. For this purpose it is useful to consider the nonlinear Hamiltonian 
of the system and to seek new terms which have positive contribution to the Hamil- 
tonian. For simplicity we restrict our consideration to terms depending on 2; = (VTr)^ 
only. 

To the leading Lagrangian Lq — {X — 1)^/2 — ^(V^Tr)^, let us add a new term of 
the form —F{z): 

1 



■{X-iy -a{V'T:y -F{z), X-1 = 2ti 



(V7r)^ (4.16) 



The Hamiltonian, or the generator of the unbroken diagonal part of (time translation 
symmetry) x (the shift symmetry) , is 

1 



H 



(4.17) 
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Fig. 13. The result of numerical simulation with a = 2.0 x lO^"', Gq = 1.0 x 10^, 
.^0 = 1-0 X 10^^ and 6 = 5.0 x 10~^. The left figure shows C,- The right figure shows 
/. It is easy to see in the left figure that |,^| larger than is cut off by the gauge 
coupling. It is also easy to see in the right figure that the implosive accumulation of 
incoming waves is supported by the repulsive electric force (/ < 0) near the center. 



The second term —z'^/2 has negative sign and, thus, the Hamiltonian would be un- 
bounded from below without new terms. Actually, this is the origin of the nonlinear 
instability. Hence, in order to stabilize the system we would like to make the Hamil- 
tonian bounded from below by choosing the new term appropriately. This is indeed 
the case if F{z) increases faster than for large z. On the other hand, F{z) must 
not change physics below cutoff M and, thus, must be turned off for small z. In 
particular, it must at least satisfy -F'(O) = 0. 

Just to show that this argument based on the Hamiltonian works, let us consider 
a simple choice 

o/ _ r /r - /9\ 1 

(4.18) 



F(z) = ^z^ 



2 



tanh (^y^j + 1 



where /3, 7 and 6 are positive constants, and F{z) = Al^{^F\z^)' . We have 
specified F{z^ instead of -^(2;) just for practical purpose: since the equation of motion 
includes the term Vtt- V2;/22;, Fi^z^ is easier to be implemented in numerical codes 
than F(^) itself. Results for the perfectly planar symmetric case with and without 
the new term are shown in Fig. 15. It is easy to see that the would-be caustics at 
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t: r 

Fig. 14. Comparison of the regular result for Gq = 1.0 x 10^ with the singular result 
for Gq = 0. The left figure shows a closer look at the regular behavior of ^ shown in 
Fig. 13. The right figure shows the singular behavior of ^ for Gq = 0. The value of 
a and the external gravitational potential $ are the same as before. Note that the 
scale of hight is different for these two figures. 

t ~ 7.5 is resolved by the new term. 

In this subsection we considered simple models of new degrees of freedom and new 
terms arising near the perfect caustics. As stated in the beginning of this subsection, 
they must be taken into account since the usual leading spatial derivative term V^tt 
vanishes for the perfect caustics and the spatial derivative expansion in this back- 
ground should start differently from that in the vr = background. We have seen 
that turning on the gauge coupling of gauged ghost condensation or /and the new 
term resolves the singular behavior. We must, however, consider the simple models 
presented here just as an existence proof of systems without the singular behavior. A 
detailed analysis of the consistent derivative expansion and/or more realistic models 
of UV completion are necessary as a future work. 

5 Bounds 

In this section, we consider the bounds on the ghost condensate. 
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Fig. 15. The evolution of X — 1 in the perfectly planar symmetric case with and 
without the new term. The left figure shows the evolution with the new term with 
13 = 0.04, 7 = 150 and 6 = 0.001. The right figure shows the evolution without the 
new term. The ghost length scale Lghost is set to be 0.01 and the initial condition at 
t = is set by vr = —0.2 x exp(— x^), X — 1 = 0. It is evident that the would-be 
caustics at t ~ 7.5 is resolved by the new term. 

5.1 Jeans Instability 

For M > 10 MeV, the Jeans instability time is shorter than the lifetime of the uni- 
verse, and we must consider the effects of this instability. We have seen that the non- 
linear effects dominate near interesting gravitational sources, but the linear dynamics 
still controls the behavior of the system for sufficiently weak ghostone amplitudes. In 
the linear regime, fluctuations with wavelength A > Lj grow on a time scale 

r~TjA, (5.1) 

where 

L,^% T,^Mk (5.2) 

are the Jeans length and time scales. Wavelengths of order Lj become unstable first, 
and longer wavelengths take longer to grow. Since fluctuations on wavelength shorter 
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than Lj are stable, we expect the minimum size of a positive or negative energy region 
to be Lj. On the other hand, the maximum size is determined by requiring that the 
time scale r above be shorter than the Hubble time. Hence, a positive or negative 
region can grow within the age of the universe if its size L is in the range 

Lj<L< L^^, (5.3) 

where 

M [ M \ , , 

^ Mpi^ ^ ^« llOOG^j • ^^-^^ 

The unstable modes grow at least until nonlinear effects become important. This 
happens for tt > tTc, where 

\2 

TTc ~ — . (5.5) 

or equivalently S > Sc with 

It is reasonable to assume that the nonlinear effects cut off the Jeans instability at 
this critical amplitude. This mechanism will fill the universe with regions of positive 
and negative ghostone field with amplitude of order iS^ and the size in the range 
Eq. (5.3). Since E is a conserved charge, there will be equal amounts of positive and 
negative E. 

The sun's Newtonian potential triggers the Jeans instability of the ghost conden- 
sate and, thus, it is expected that there be a positive or negative region around the 
sun. This is justified if the 'aether' is efficiently dragged by the sun and we now argue 
that this is indeed the case. To do this, it is useful to work in the rest frame of the 
sun. Far from the sun, the aether is moving with constant velocity v ~ 10^^, but near 
the sun the velocity field will be distorted by the presence of the sun. By using the 
fluid picture of the ghostone fleld, we estimate the effect on a ffuid particle with speed 
V and impact parameter r. The fluid particle will be a distance of order r away for a 
time At ~ r/v, so the change in the particle velocity in the impulse approximation is 

Rs r Rs , . 

At! ~ — • - ~ — , (5.7) 
r V vr 

where i?s is the Schwarzschild radius of the source. Thus, the change in the velocity of 
a fluid particle becomes comparable to or greater than the initial velocity if r < rdrag, 
where 

rdrag ~ (5.8) 
V 
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For our sun, rdrag ~ lOi?©, so the dragged region extends outside the solar radius.^ 

We require that the absolute value of the mass of the lump with the critical density 
Pc and the size I/max be at worst less than the solar mass: 

PcLl^ ^ Me- (5.9) 

This requirement gives the bound 

M ^ 10^ GeV. (5.10) 

Since the high power of M (the l.h.s. oc M^) is involved in Eq. (5.9), a more stringent 
requirement on the mass of the lump will not substantially improve the bound. 



5.2 Twinkling from Lensing 

We have argued that if M > 10 MeV, then the Jeans instability fills the universe with 
regions of positive and negative energy of size L > Lj ~ M-pi/M"^ with energy density 
Pc ^ M^/Mp^. This will happen everywhere, in particular in the voids between 
galaxies. Any light that travels to us from far away will therefore be lensed by these 
positive and negative regions. These positive and negative energy regions move, 
because the local rest frame of the lensing regions is different from that of our galaxy, 
so the result is that the observed luminosity of any point source will change with 
time. This is similar to the twinkling of the stars in the night sky caused by time 
dependent temperature differences in the atmosphere. In this subsection, we work 
out the bounds on the ghost condensate from this effect. 

Suppose that the universe is filled with regions of positive and negative energy 
with size L and density Pc- A light ray traveling through such a region will lens by 
an angle 

~ $ ~ ^ ~ (5.11) 

If a light ray travels a distance (i 3> L, then it will undergo N d/L uncorrelated 
lensing events, so the total angular deviation will be enhanced by a N^/"^ random walk 
factor: 



We see that the largest angular deviation comes from the largest L and largest d 



^This radius is still much less than the orbital radius of Mercury. 
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The size of L is limited by the time for the Jeans instabihty to form as in (5.3). 
If the source is the cosmic microwave background, then d ~ Hq^ and we obtain 

for the largest regions with the size L ~ L^nax- The high power of M makes the 
precise experimental limit on A^cmb irrelevant, and we obtain the bound 

M < 100 GeV. (5.14) 

For M ~ 100 GeV, the size of the largest critical region is L ~ 10^^ cm, approxi- 
mately the radius of the sun. The local velocity of these regions relative to our galaxy 
is of order 10"^, so the time scale for one of these regions to cross the line of sight is 
of order a day, which is therefore the time scale of the variation. 

If there is a distant astrophysical source that is observed to shine with very little 
time variation, it may give a competitive bound. But given the high power of M 
involved, it seems difficult to improve on this bound significantly. 

5.3 Supernova time- delay 

Gravitational lensing considered in the previous subsection induces a time-delay for 
light-rays coming from far distances. With this time-delay effect, observed supernovae 
should be older than they appear. Thus, this effect would change the estimate of dark 
energy by observation of Type la supernovae. Since the determination of the dark 
energy by supernovae observation is known to be consistent with the WMAP data, 
we require that the time-delay is sufficiently shorter than the total time: 

Note that the precise experimental limit on the At/tis irrelevant because of the higher 
power of M involved in the l.h.s. Prom this we obtain the bound 

M ^ 10^ GeV. (5.16) 

5.4 Nonlinear Energy Loss 

A gravitating object moving relative to the preferred frame will lose energy. As dis- 
cussed above, the ghostone field outside a gravitating object is in the nonlinear regime 
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for essentially all gravitationally interesting objects and interesting values of M. In 
this subsection, we estimate this energy loss in the nonlinear regime. Interestingly, 
we obtain parametrically the same result as the linear calculation, and so this is a 
small effect in the real world. 

It is simplest to work in a frame where the ghost fluid is at rest. As the gravitating 
object moves through the medium, its gravitational potential excites the surround- 
ing ghost fluid, and therefore the object loses energy. This picture is very similar 
to the energy loss of charged particles in a medium, treated in standard texts on 
electrodynamics. In fact, the particle picture makes the analogy very close. 

Consider then a gravitating object characterized by its Schwarzschild radius Rs = 
M^/Mpi moving with speed v in the +z direction relative to the preferred frame. 
We estimate the effect on a particle with impact parameter r using the impulse 
approximation in order to see the response of the ghost condensate in the particle 
picture. The moving object will be a distance or order r away for a time At ~ r/v, 
so the change in the particle velocity in the impulse approximation is 

~ . ^ ~ (5.17) 

r V vr 

In terms of the ghostonc field, this means that we generate a disturbance with Vvr ~ 
Av. To estimate the energy density associated with this, we must remember that 
there is no modification of gravity if we neglect the Q;(V^7r)^ term. The change in the 
ghostone field energy is therefore 

Ap ~ aM^V'n)^ ~ (5.18) 
The total energy transferred to the ghostone fiuid is therefore 

'Jdrr-^. (5.19) 



^ dE f, aM^Ri 
h ~ V—— V Idrr- 
dz 



This integral is dominated at small r, so we have 

E^'^. (5.20) 

vr ■ 

^' mm 

We now consider the short distance cutoff r-^am- One cutoff comes from the use 
of the impulse approximation. This is valid only if the particles moves a distance 
smaller than r in the time At: 

<^tr-^i^y^r, (5.21) 
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which gives 

Rs 



' mm o 
V 



(5.22) 



For r ~ rmin we have Av ~ v, so the maximum velocity of particles with r < Tmin is 
Aw ~ V. This shows that the contribution of particles with r < rmin is suppressed, 
justifying the use of r^^i^ as a cutoff. 

Another cutoff is provided by the finite size R of the object itself. This is smaller 
than Tmin given above provided that 

§ Z v\ (5.23) 

The left-hand side is the surface gravity of the object, and the right-hand side is of 
order 10^^ for a typical astrophysical object. For example, the sun has a surface 
gravity of order 10"^, while the earth has a surface gravity of order 10~^. 

Putting this together, we see that for an object of sufficiently strong surface gravity 
(Eq. (5.23)), the energy loss is 

E ~ aM'^v^. (5.24) 

Note that the dependence on the mass of the object (as well as the Planck mass) 
has dropped out. Eq. (5.24) is the same parametric formula obtained from a linear 
calculation, and completely negligible for macroscopic objects. If the surface gravity 
is weak, the energy loss is even smaller. 



6 Ghost Dark Matter 

As shown in Ref. [5] and reviewed above, the expansion of the universe drives E — > 0. 
If E is nonzero early in the history of the universe, then the ghostone energy density 
evolves with the scale factor a as 

p ~ E ~ a-^ (6.1) 

with p ~ p^/M^ <^ p for p <^ M^. (Since the energy density of the universe at 
matter-radiation equality is ~ (1 eV)'', this requres M > 1 cV.) Thus, the ghostone 
energy and pressure redshifts just like cold dark matter. Note that since E is a 
conserved charge density, a nonzero value of E early in the universe corresponds to a 
nonzero charge asymmetry, similar to the baryon asymmetry in the universe. In such 
a scenario, the amount of dark matter is determined by very different microphysics 
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from the conventional freeze-out of particle dark matter. In particular, if we want 
to generate a nonzero density of E charge from an initial condition with no charge 
asymmetry ('ghostogenesis'), the theory must contain interactions that violate the 
shift symmetry. We will not discuss the required initial conditions for ghost dark 
matter further, but instead turn to the viability of this scenario. 

As already noted in Ref. [5] , small fluctuations of the ghostone field with 5p/ p <^ 1 
also behave just like cold dark matter. This can be seen from the fact that the equation 
of motion in both cases follows from the conservation of the stress-energy tensor for 
linearized scalar fluctuations, which is the same for the ghostone field as for cold dark 
matter. Using the results of the present paper, wc can sec explicitly that the equations 
of motion for the ghostone mode are equivalent to a perfect fluid with equation of 
state p = /2M^ even in the nonlinear regime bpj p 1. 

However, the ghostone fluid behaves very differently from cold dark matter in 
the nonhnear regime. To see this, let us recall the basic facts about the nonlinear 
dynamics of cold dark matter. An overdense region of cold dark matter starts gravi- 
tational collapse, and decouples from the universe expansion. The cold dark matter 
particles do not reach the center of the overdense region because they carry angular 
momentum, and there is no efficient mechanism for them to lose it. The trajectories 
of cold dark matter particles gcncrically cross, but nothing special happens because 
the particles have negligible interactions. The over-dense region therefore virializcs, 
becoming effectively a cloud of dark matter particles orbiting in the gravitational 
potential produced by all the other dark matter particles. At any given point in the 
overdense region, there are streams of dark matter particles orbiting in many differ- 
ent directions. Numerical simulations show that this kind of dark matter is quite 
successful in reproducing many of the observed features of structure in our universe, 
although there are open questions and possible discrepancies. (For a recent review, 
see e.g. Ref. [16].) 

In the case of ghost dark matter, the trajectories of the 'particles' in the fiuid 
picture cannot cross, since the theory is microscopically a scalar field. The pressure 
term will provide a repulsive force in the over-dense regions, and a simple argument 
suggests that it is more important than higher derivative terms such as the Q;(V^7r)^ 
term in the action. In a region where ghostone energy is the dominant gravitational 
source, the gravitational potential is given by 



(6.2) 



where 



M2 



1 



(6.3) 



m = 
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The equation of motion of the fluid particles is 

^=-V($ + E), (6.4) 

— * 

where the VE term gives rise to the pressure force. Prom this we see that the pressure 
force will dominate over the gravitational forces if S > $. In the case where the 
dominant source of gravity is E (as in a halo of ghost dark matter), Eq. (6.2) shows 
that VS > V$ for lumps smaller than ~ Lj. This means that the pressure force 
will act to smooth out any inhomogencitics on scales smaller than Lj. 

To see that the pressure force can indeed balance gravitational forces, we can 
consider spherically symmetric static configurations, i.e. 'ghost stars.' The force- free 
condition is V($ + S) = 0, which implies E = — $ + constant. If the constant is 
nonzero, then $ ~ at large r, which can be discarded because it has growing 
energy density at infinity. We therefore have E = — and Eq. (6.2) gives a solution 

E = -$ = Eo^^, (6.5) 
mr 

where So = E(r = 0). For this solution the vr's gradient grows, but its contribution 
to the energy density is suppressed by M^/Mpi and thus negligible. Actually, the 
gradient of n at time t after the formation of the 'ghost star' is estimated as 

V^TT ~ m'^nt ~ m^Eot ~ ^^^^^ > (6-6) 

Mpi tdyn 

where tdyn ~ ^Gnpq dynamical time for the 'ghost star' and po — M^TIq is the 

central density. Thus, the energy due to the gradient is estimated as 

Ap~aM2(VV~a^Po(— J • (6.7) 



Hence, Ap/po remains suppressed by the factor M^/Mpi and is negligible for t ~ tdyn- 
This means that the 'ghost star' can be considered as a quasi-static configuration. 
This solution has alternating regions of positive and negative S, and therefore does 
not describe a normal star. However, the region r < vr/m has positive S (if Eq > 0) 
and describes a region where gravity and pressure forces balance. In fact, we can 
easily check that the solution in this region is stable. The Lagrangian for fiuctuations 
about this solution is 



1-2 ly. 



(6.8) 
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The equation of motion for spherical perturbations is 



^0 

TT = -dr [r sin(mr)9r7r] . (6-9) 



Separating variables as tt = ■7f(r)e multiplying the equation of motion by r^Ti", 
integrating from r — Otor — Tr/m and using integration-by-parts, we obtain 



u;^ / ['^(^)] r^dr = — / rsm(mr)dr. (6.10) 

Jo m Jo 

This implies that o;^ > unless n' — everywhere in the interval of integration. 
Hence, the positive E region near the origin is stable. Although this solution is 
not realistic, it gives an explicit example where pressure can cancel gravitational 
attraction. 

We see that in regions where E > is the dominant source of gravity, the pressure 
force tends to smooth out inhomogeneities on scales smaller than Lj. Let us ask 
whether the Q;(V^7r)^ term is important in these regions. In the solution without the 
Q;(V^7r)^ term, we have 

(6.11) 

Prom this, we can see that the Q;(V^7r)^ term is neghgible for 

L > (6.12) 
This is satisfied at L ~ Lj provided that 

E Z (6.13) 

Note that this is the critical amphtude Ec discussed in subsection 5.1. If we assume 
that the ghostone energy density p ~ M^E accounts for all dark matter in the uni- 
verse, then Eq. (6.13) is always satisfied in the overdense regions for M ^ 10 MeV. 
In that case, it is very plausible that the pressure is sufficient to smooth out the 
would-be caustics that occur during structure formation in our universe. 

We emphasize that we arc not claiming that the pressure prevents caustic singu- 
larities from forming during structure formation. For example, a collapsing E > 
region may 'bounce' and leave behind a region of E ~ (a void) or even E < 0. In 
these regions, the singularities discussed in the previous section may become impor- 
tant. Also, a very over-dense region may lead to caustics in a less dense region nearby. 
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However, even if the a(V^7r)^ term becomes important in some small regions of the 
halo, leading to negative energy regions, we argued in the previous section that the 
negative energy regions tend to shrink and carry negligible energy and they happen 
at very small scales. Therefore, we do not expect them to be important in the overall 
distribution of energy. 

It would certainly be very interesting to perform numerical simulations of the non- 
linear dynamics of self-gravitating ghostone halos to see whether ghost dark matter 
can explain the structure of the universe we observe today. 

7 Conclusions 

Infrared modifications of gravity have potential connections to the puzzles wc arc 
facing in cosmology such as dark energy and dark matter. Such theories can be viewed 
as gravity in the Higgs phase. The simplest possibility with a single extra degree of 
freedom which nonlinearly realizes the broken time diffeomorphism is described by 
ghost condensation — a scalar field with a time dependent vacuum expectation value, 
(f) (X t. In a previous paper [5], we have shown that a healthy low-energy effective 
theory can be constructed following the symmetry requirement. It is stable for small 
fluctuations in the absence of gravity (Mpi — > oo). When gravity is turned on, there 
is an infrared instability at long wavelengths and large time scale similar to the usual 
Jeans instability for ordinary matter in the presence of gravity. This potentially 
provides the strongest constraint on the UV cutoff of the effective theory. 

The leading nonlinear interactions in the effective theory of ghost condensation 
are also fixed by the symmetry. By power counting, they are irrelevant operators so 
they are unimportant for small fiuctuations about the ground state. However, they 
can become important in the presence of large classical sources. In this paper, we 
show that for realistic astrophysical gravity sources, the time scale for the nonlinear 
terms to become large is shorter than the time scale of linear propagations. Therefore, 
the nonlinear effects are dominant near these astrophysical bodies and this paper is 
devoted to studying these effects. 

Wc show that there is a simple picture for the nonlinear evolutions of the 'ghost 
field.' They arc simply described by the geodesies in the gravitational field of the 
source. The geodesic picture necessarily breaks down when geodesies cross and caus- 
tics form. At these points the higher derivative terms become important. Prom the 
numerical simulations we see that the generic caustics can be resolved by the leading 
higher derivative terms in the linear lagrangian. However, there exist perfect caustics 
which cannot be stopped by the leading higher derivative terms. Their resolution can 
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only be addressed by the physics above the UV cutoff, beyond the regime of vahdity 
of the low-energy effective theory. A potential worry is whether some catastrophe 
may be caused by the negative energy region formed due to the caustics. We show 
that both analytically and in numerical simulations that the negative energy region 
shrinks with the total amount of energy inside the region staying approximately con- 
stant, so it will not swallow the whole space. In addition, for small enough Af the 
total energy inside each of these region is extremely small by astrophysics standards 
so it is not expected to give large anomalous gravitational effects. As long as such 
regions can be resolved in a full theory, we do not expect disastrous effects or strong 
bounds coming from these caustics. 

In fact, the nonlinear effects may even relax the bound on the UV cutoff of the 
theory. In the linear approximation, we obtained a bound of M ;^ 10 MeV from 
requiring the time scale of the Jeans instabihty to be longer than the age of the 
universe. The nonlinear effects can cut off the growth of the Jeans instability and 
hence removes the above bound. The total energy in a lump can only grow to ~ 
Mpi which is minuscule before the nonlinear dynamics cut off the growth. Now the 
question is how the positive and negative energy regions change the standard picture 
of our universe. We found that light rays coming from far distances lens by those 
lumps. By requiring that the random walk of light rays due to the sequence of lensing 
events docs not smear out the observed CMB anisotropics, we found a slightly weaker 
bound of M < 100 GeV. 

Nonlinear dynamics will also play an important role if we try to imagine that 
dark matter is made of ghost condensate. We have shown that a coherent deviation 
of ghost condensation from the ground state redshifts like matter, and in the linear 
regime the growth of the perturbations is identical to the usual cold dark matter. 
Whether it can give rise to the correct structure formation in the universe depends 
on the nonlinear evolution. This is an interesting but complicated question, which 
may require detailed simulations. It is left for future studies. 

In this paper, we have begun the investigation of the nonlinear dynamics of ghost 
condensation as a Higgs phase of gravity, which exhibits a rich structure. While 
we have obtained a qualitative picture of some of the dynamics, a more complete 
understanding is clearly desirable and deserves further investigations, as required for 
making more concrete connections with various cosmological questions. 
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Appendix: Higgs Phase of Full General Relativity 



A.l Fluid Picture 

We now show that the fluid picture derived in the main text holds in full general 
relativity, i.e. allowing relativistic motion and strong gravitational fields. In order to 
remain within the realm of validity of the effective theory, we must restrict ourselves 
to small ghostone ampltiudes, i.e. 

S<L (A.l) 

Neglecting the a term in the action, the equations of motion of the ghostone field can 
be written 

V^J^ = 0, (A.2) 

where 

J" = Eu", (A.3) 

and 

u^" = V^0. (A.4) 

We cannot interpret as a 4-velocity because = X ^ 1 in general. However, this 
effect is C(S^), as can be seen by defining a normalized 4- vector by 



1 + 0(E^)1. (A.5) 
We can therefore interpret 

u^" = (0, -V(/.) (A.6) 
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as the 4-velocity of the fluid, up to corrections. 

It now remains only to satisfy the relation between S and 0. Differentiating the 
relation 

^''"V^0V,0 = X = 1 + 2E, (A.7) 

we obtain 

V^E = iV^(V''0V,0) = W<PV^V,<P = Wcl>V,V^ci> (A.8) 

or 

u'VpU^ = V^S. (A.9) 

This has the form of the geodesic equation of a particle subject to an additional 
gradient force with potential S. For the special case E = 0, the fluid particles follow 
geodesies. The velocity fleld satisfles the curl-free condition 

V[^u,] = 0, (A.IO) 

implying that the flow is irrotational. 

Conversely, any family of particle trajectories satisfying Eq. (A.9) with initial 
conditions 

u^l = given, V[^ii^]| = (A.ll) 

on some initial value hypersurface defines a solution of the ghostone equations of 
motion. To see this, note that the vanishing of V[^Mj^] is preserved by the equation of 
motion Eq. (A.9). This can be seen from the identity 

u'Vp{V[^u,]) = (V[.Mp])VX - (V[;.Mp])V.M^ (A.12) 
which is derived as follows: 

= -VXV[pH,] - VXV.Hp + V^V.E - Rp^.^^uPu". (A.13) 

Antisymmetrizing in ji and v then gives Eq. (A.12). Therefore, the velocity field 
can be written as a gradient Ufj_ — V ix<t>- The field now obeys the equation of motion 
Eq. (A.2), with 

51^'' = X + constant. (A.14) 
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The integration constant can be seen to vanish by imposing appropriate boundary 
conditions at infinity. 

This means that many of the results derived in the main text hold in full general 
relativity. In particular, in regions where the solution is given by E = 0, the local 
preferred frame given by is a freely falling frame. In this sense, the aether is 
'dragged' even by strong gravitational fields. 

A. 2 Negative Energy Region Shrinks 

We now show that any S < region shrinks in full general relativity, neglecting the 
a(V^7r)^ term in the equation of motion. 

Now let us consider a hypersurface E = defining the boundary between a nega- 
tive energy region (E < 0) and the surrounding positive energy region (E > 0). On 
the boundary E = 0, the equation of motion is reduced to 

a^EM^ = 0. (A. 15) 

Now let us calculate the acceleration a'^ = u'^'VpU'^. By differentiating the definition 
of E, we obtain 

9^E = uPV^Up = u^VpU^, (A. 16) 

where we have used V[juMp] — 0, and thus 

= 9^E. (A. 17) 

With this expression it is easy to know the direction to which the boundary hyper- 
surface evolves: 

a'^V^E = 9^E9^E < 0. (A. 18) 

Namely, the boundary hypersurface evolves towards the region with negative E. In 
other words, the negative energy region shrinks. 

A. 3 Relativistic Caustic 

We now show how the caustic solutions discussed in subsection 3.5 can be found in an 
elegant way in the relativistic formulation. We consider the case of 1 + 1 dimensions, 
which is relevant for the generic caustic. We have seen that the spacetime surface 
(f){x) = is perpendicular to the family of geodesies. We can therefore use Gaussian 
normal coordinates based on this surface. In this coordinate system, (p is the 'time' 
variable, and the metric in the {x, t) subspace is 

ds'^ = #2 - A\(P, a)da^, (A. 19) 
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where the spatial variable a parameterizes the position along the = constant sur- 
faces. Without loss of generality, we can choose ^4(0, a) = constant, since this corre- 
sponds to choosing the 'space' coordinate on the surface = to be proper distance. 
The geodesies perpendicular to the surfaces of constant are simply a — constant in 
these coordinates. 

The condition that the metric Eq. (A. 19) is flat is 



dlA = 0. 



(A.20) 



(For example, one can check that this guarantees the vanishing of the Riemann ten- 
sor.) Imposing ^4(0, cr) = 1, the general solution is 



A{cj>,a) = [l 



(A.21) 



where 0c(c) is a function that specifies the distribution of velocities on the surface 
= 0. Note that the coordinate system breaks down when = 0c(c); this is a sign 
of the caustic singularity. For a caustic singularity 

In this coordinate system, it is trivial to compute invariants such as 



'-9 



-A. 



\d,A 



1 

T 



(A.22) 
(A.23) 



where 



r = - 0c 



is the time to the caustic. Similarly, 
1 



(A.24) 



□ ^0 



-9 



-g g'^'d, □ 0] = [Ad^ 00]- [A-'d, □ 0] (A.25) 



0c M.r + 0c0'c'] 



+ 



Expanding the function 0c(cr) for small a, we obtain 

1 mm 



(A.26) 



(A.27) 



This agrees precisely with the results of subsection 3.5 in the nonrelativistic limit, 
with the identifications a — Xq, (l)c{cr) — tc, and t — t — T. 
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